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Abstract. The present paper describes the effect of magnetic pick-up coil transfer 
functions on mode number analysis in magnetically confined fusion plasmas. Magnetic 
probes mounted inside the vacuum chamber are widely used to characterize the mode 
structure of magnetohydrodynamic modes, as, due to their relative simplicity and 
compact nature, several coils can be distributed over the vessel. Phase differences 
between the transfer functions of different magnetic pick-up coils lead to systematic 
errors in time- and frequency resolved mode number analysis. This paper presents the 
first in-situ, end-to-end calibration of a magnetic pick-up coil system which was carried 
out by using an in-vessel driving coil on ASDEX Upgrade. The effect of the phase 
differences in the pick-up coil transfer functions is most significant in the 50-250 kHz 
frequency range, where the relative phase shift between the different probes can be up 
to 1 radian (~ 60°). By applying a correction based on the transfer functions we found 
smaller residuals of mode number fitting in the considered discharges. In most cases 
an order of magnitude improvement was observed in the residuals of the mode number 
fits, which could open the way to investigate weaker electromagnetic oscillations with 
even high mode numbers. 
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1. Introduction 

Magnetohydrodynamic (MHD) instabilities have a strong influence on the plasma 
performance, therefore their thorough understanding is essential for stable tokamak 

f See http://www.euro-fusionscipub.org/mstl 
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operation. Early steps of identification and analysis inclnde determining the time- 
freqnency evolntion, spatial strnctnre (characterized by mode numbers) and whether 
the perturbations are global or localized. Several basic measurements of MHD modes 
rely on magnetic pick-up coils mounted on the inner side of the vacuum chamber. 
Depending on the alignment these probes pick up different components of the magnetic 
held huctuation. On ASDEX Upgrade (AUG), Mirnov coils are placed on the vacuum 
vessel to measure the poloidal component of the magnetic held huctuations and so- 
called ballooning coils are placed closer to the plasma on the low held side to measure 
the radial component Due to their relative simplicity and compact size, several 
coils can be distributed over diherent locations within the plasma vessel, making them 
invaluable candidates for the analysis of the mode structure. 

Several methods are available to determine mode numbers of global MHD modes. 
One example is the application of Singular Value Decomposition (SVD) [1, 2], which 
decomposes huctuation data by choosing the basis vectors in the directions of maximum 
coherence. With spatially distributed measurements one can obtain both the coherent 
spatial structure and the time evolution of a mode. Another promising way of 
interpreting spatial and temporal data is the use of the Lomb periodogram, which was 
proposed by Zegenhagen and co-authors [3]. Our preferred approach, which is applied 
throughout this paper, is to apply linear continuous time-frequency transformations on 
the data and calculate the mode structure based on the phase diherences between probe 
pairs [4, 5]. A careful treatment of mode structure analysis takes into account the 
screening currents induced in the vessel wall and other in-vessel conducting components 
by a simulation of the perturbed magnetic held from an assumed current distribution 
on resonant surfaces [6, 7]. 

Differences in the phase behaviour of the transfer functions corresponding to 
different coils can introduce systematic errors in the mode number htting algorithms. In 
order to correctly characterize mode numbers, each of the mentioned methods need to 
be evaluated on calibrated signals of the magnetic probes, which takes into account the 
effect of the magnetic coil response. Systematic errors come from various sources such 
as eddy currents, uncertainties in the probe positioning and alignment, resonances in 
the probe or the measurement circuit, due to the anti-aliasing analogue hlter or even the 
analogue-digital converter (ADC). We discovered systematic errors in the mode number 
analysis of strong MHD modes in the 50 — 250 kHz frequency range on ASDEX Upgrade. 
Several methods have been previously applied to model or measure the transfer function 
of the magnetic coils, however no in-situ, end-to-end calibration of the full diagnostic 
system using an in-vessel driving coil has been done before. A remote calibration method 
has been developed on JET [8, 9]. This method determines the transfer function for the 
ampliher & digitizer section by directly injecting a known signal. The probe & cable 
section of the system was modeled as a lumped-circuit and its transfer function was 
extracted from the frequency-dependent impedance measurement performed without 
entering the torus. On TCV the response of the magnetic probes and the amplifying 
chain were measured separately before installation and a transfer function with two cut- 
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off frequencies was fitted on the product of the two evaluated transfer functions [10, 11]. 
A Helmholtz coil was used to measure the transfer function of the coils outside the 
vacuum chamber on H-INF heliac [12], On TFTR, only the electronics of the Minov 
coils were calibrated [13]. 

This paper presents the evaluation of the first in-situ, end-to-end calibration of a 
magnetic pick-up coil system which was carried out using an in-vessel driving coil [14]. A 
phase correction function was derived from the coil response measurements and it was 
applied to improve mode number determination on AUG. We found that the effect 
of the phase differences in the pick-up coil transfer functions is most significant in 
the 50-250 kHz frequency range, where the relative phase shift between the different 
probes can be up to 1 radian (~ 60°). The implications of this systematic error 
and the corresponding correction is demonstrated on strong Toroidicity-induced Alfven 
Eigenmodes (TAEs) observed in AUG discharges, where a significant improvement in 
toroidal mode number fitting is achieved. 

The paper is organized as follows. In section 2 we present the mathematical 
background of our mode number analysis and the set of magnetic pick-up coils used. 
Systematic errors observed in magnetic measurements are described in section 3. 
Section 4 details the in-situ transfer function measurements on ASDEX Upgrade. The 
implementation of the correction, including examples, is described in section 5; followed 
by the conclusions in section 6. 

2. Analysis principles and measnrement set-np 

Our mode number analysis uses linear continuous time-frequency transforms [15]. These 
are calculated by expanding the signal f{t) on the basis of families of so-called time- 
frequency atoms: 

+ 00 

Tf{t,uj) = {f,gt,u;) = j f{t')gl^{t')([t' , (1) 

— CO 

where gt,uj is a time-frequency atom, whose energy is well localized in both time and 
frequency. Variables t and u are the time and frequency indices of the atom identifying 
its position on the time-frequency plane, and the * represents the complex conjugation. 
The energy density distribution can then be calculated by taking the absolute value 
squared [15]: 

E/(^,n;) = |T/(^,a;)^ (2) 

Our mode number calculation process is based on the phase of the cross-transform: 

= arg{T/,(t,a;)T/;(t,a;)} , (3) 

where fk and fi represent the signals of probes placed in different positions. 

In this paper, we demonstrate the mode number determination method by using 
Short Time Fourier Transform (STFT) [15]. Note that this method also works with 
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continuous analytical wavelet transform [15]. STFT is a continuous time-frequency 
transform where the family of time-frequency atoms are generated by shifting a real and 
symmetric window g{t') in time {t) and frequency (u;): 

gtAt')=e^^^'g{t'-t) , (4) 

which gives a constant time-frequency resolution on the time-frequency plane. Following 
from (1) STFT can be written as 

Sfit,u) = {f,gt,co) = 

+0O 

= j f{t')e~'‘^^'g{t' - t)dR , (5) 

— OO 

where g{t') is a Gaussian function in our case. The energy density distribution dehned in 
(2) calculated from STFT is called a spectrogram, which is ideal to investigate the time- 
frequency evolution of non-stationary modes [16, 17]. In order to apply a continuous 
transform on discrete time signals, the transform has to be discretized in a way to avoid 
the degradation of the time-shift invariance property of the transform [18]. 

Hereinafter, we focus on toroidal mode number analysis on AUG, derived from 
the measurements of six ballooning coils placed close to the separatrix (~ 10 cm) at 
the outer midplane in different toroidal, but almost identical poloidal positions. These 
ballooning coils are rectangular spirals (etched out vacuum-compliant printed circuit 
boards) which measure the changes in the radial component of the magnetic held. The 
vanishing dimension in the direction of the held makes them extremely sensitive for 
short wavelength perturbations. The toroidal positioning of these probes is illustrated 
in hgure 1. The measurement set-up on AUG contains ADGs which have a 2 MHz 
sampling frequency and they include a 3dB low pass Bessel hlter with 512 kHz cut- 
oh [19]. 

The measurement of magnetic held huctuations produced by MHD modes allows 
the reconstruction of their mode structure. Harmonics of global MHD eigenmodes are 
generally assumed in the following form [3] : 

Am,n(P,6'*,0,t) = 

= Am,n{.P, 0*) exp{i(m6'* + ncf - tot)} , (6) 

where A is the observable quantity (in our case the radial magnetic held perturbation 
B-f) as a function of time t and the so-called straight held line coordinates (p, 9*, 0) [6], 
A{p,6*) is the radial eigenfuntion and u is the mode frequency. The mode structure 
of the harmonic is characterized by the m poloidal and n toroidal mode numbers. As 
follows from (6), in the case of a single, pure sinusoidal, globally coherent mode at 
hxed {p,9*,t), but at diherent toroidal angles, the relative phase between Ara,n{.4>k) and 
Am,n{4>i) is proportional to the relative toroidal angle (j)k — 4>i and the ratio is the toroidal 
mode number n\ 

Vki\ = arg{A^,„(0fc)} - aig {A^^n{(j)i)} = 

= n{(j)k - 4>i) ■ 


( 7 ) 
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Figure 1. A top-down view of ASDEX Upgrade indicating the position of the 
ballooning coils in the toroidal array. The arrow marks the B31-14 probe which is 
used as a reference. 

Thus, in order to characterize the spatial structure of the mode, one has to determine its 
phase in different locations {4>i). Equation (7) shows that the slope of the linear curve, 
htted on relative phases between all pairs of signals as a function of the relative probe 
position, gives the toroidal mode number n. Our method searches for the integer mode 
number value where the residual Q of the fit is minimal: 

Q = ~ ) ( 8 ) 

V 

where ip-p is the phase between pairs of signals, (fp is the relative probe position, n is the 
toroidal mode number, the sum is executed on the V signal pairs and ||.|| is the norm 
obtained by taking the optimum shift of (pp by 2Trz, where ^ is an integer number. The 
sign of the mode number determines the direction of propagation in the device frame. 
If the helicity (which is determined by the direction of the toroidal magnetic field and 
the plasma current) is known, the sign of the mode number can be related to the ion or 
electron diamagnetic drift direction. 

A time-frequency resolved map of the mode numbers can be generated by 
performing the linear fit in each time-frequency point [4]. Note that this result is 
only valid, where global modes exist. The time-frequency based mode number analysis 
presented above is implemented in NTI Wavelet Tools custom data processing toolbox§. 

In principle, the calculated cross-phase in (3) is a reasonable approximation of the 
relative phase of the mode, however, if the phase information carried by the measured 
signals suffer from systematic errors, the mode number analysis can yield erroneous 
results. 

§ Developed and maintained in the Institute of Nuclear Techniques (NTI), Budapest University of 
Technology and Economics, https://deep.reak.bme.hu/projects/wavelet 
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3. Systematic errors in magnetic measurements 

The systematic errors can be easily recognized in the presence of strong, coherent MHD 
modes. In such a case, one would expect that the calculated relative phases (eq. (3)) 
between all pairs of signals would lie on a straight line as a function of the relative probe 
position. As we demonstrate in the following example, systematic deviations from the 
best fitting line were observed. In our example we follow the time evolution of a strong 
TAE illustrated in hgure 2. Since 6 probes are available in the toroidal array, 15 pairs 
of coils can be selected. 


kHz 

kHz 

kHz 

kHz 

kHz 

kHz 


1.25 1.30 1.35 1.40 

Time [s] 

Figure 2. Magnetic spectrogram of strong coherent TAE harmonics with different 
toroidal mode numbers. The frequency of the mode with n = —4 toroidal mode 
number (as calculated in figure 3) is indicated with black dashed line. The time- 
frequency points indicated with colours were chosen for mode number analysis. 

We trace the time evolution of the mode frequency using a ridge-following 
algorithm [20, 21], the result of which is called the ridge of the STFT transform. The 
frequency ridge of the mode with toroidal mode number n = —4 is indicated with 
dashed black line on figure 2. We mark the time-frequency points selected for mode 
number analysis by the cross section of the time points indicated with white vertical 
lines on hgure 2 and the ridge of the mode. The relative phases calculated between 
all signal pairs are plotted as a function of relative probe position on hgure 3a. The 
diherent colours represent the time points corresponding to hgure 2. The linear ht on 
the relative phases is performed in each time instant trying to ht integer mode numbers 
from the interval [—100,-1-100]. The residuals (dehned in (8)) of the ht at 1.32 sec are 
shown as a function of the mode number in hgure 3b. In the interval [—20, -|-20] the 
best htting toroidal mode number is n = —4, which also holds in the other investigated 
time points. However, we observe a systematic deviation of the relative phases from the 
n = —4 line (plotted with black solid lines on hgure 3a). The deviation of the relative 
phases from the n = —4 for independent channel pairs as a function of time is shown in 
hgure 4. Such systematic deviations were observed in all investigated cases. 
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Figure 3. (a) The relative phases between all pairs of signals are plotted as a function 
of relative probe position. The different colours present the different time points 
depicted in figure 2. The best fitting line with n = —4 slope is plotted with black 
solid lines. For the helicity of this discharge, the negative mode number means that 
the mode propagates in the ion diamagnetic drift direction, (b) The residuals (defined 
in (8)), which quantify the differences between the fit and the measurements, are shown 
as a function of the mode number. In the [—20, +20] range n = —4 is the best fitting 
value. 
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Figure 4. The deviation of the relative phases from the n = —4 for independent 
channel pairs as a function of time. The systematic behaviour of the deviation is 
clearly visible. 


Such systematic errors are likely to arise from the differences in the phase behaviour 
of the transfer functions corresponding to different magnetic coils. In the next section 
we present the results and implications of in-situ magnetic coil response measurements 
performed on ASDEX Upgrade. 















Reducing systematic errors in time-frequency resolved mode number analysis 

4. Coil response measurement on ASDEX Upgrade 


We modeled the measurement circuit of a magnetic probe as a linear single-input single¬ 
output (SISO) system. The input of the system is in principle the time derivative of 
the radial magnetic held in the centre of a ballooning coil and the output is the signal 
registered by the ADC. The frequency-dependent transfer function of a linear system at 
a given frequency uj is 


H{uj) = 


(M 

I{u) 


(9) 


where 0{u) and I{u) are the Fourier transforms of the output and input signals. 

In order to determine the transfer function of the ballooning coils an in-situ 
measurement was performed on ASDEX Upgrade [14]. Since this measurement was 
carried out inside the fully equipped tokamak device on the installed magnetic probes, 
the result contains the effect of the coil response, the analogue hlter (a low pass, 3dB 
Bessel hlter with 512 kHz cut-oh), the ADC (2 MHz sampling frequency), the signal 
transmission cables and to some extent even the ehect of eddy currents in the in-vessel 
structures. 

The measurements from which the transfer functions of the magnetic coils are 
estimated were carried out by exciting the magnetic probes one by one with a driving 
coil placed in front of the probes. The geometry of the measurement set-up is shown in 
hgure 5. The input voltage of the driving coil came from a constant amplitude function 



Figure 5. The geometry of the measurement set-up from which the transfer functions 
of the magnetic coils are estimated. The red pick-up coil was excited by the 
green driving coil. The black supporting structure provided the consistency of the 
measurement geometry from probe to probe. 

generator, that swept the 1 kHz - 1 MHz frequency range in 4 kHz steps with a 9 second 
repetition time. The current flowing in the driving coil was measured with a current 
clamp (clamp-on AC current meter). The magnetic held at the probe is directly related 
to the current (Xd) howing in the driving coil: 

B CdXd , 


( 10 ) 
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where Cd is a constant which gives the strength of the magnetic field per nnit cnrrent and 
has a dimension T/A. Its valne is nneqnivocally determined by the geometric strnctnre of 
the measnrement set-np. Since the magnetic probe detects the changes of the magnetic 
field, we need the time derivative of (10), which corresponds to a mnltiplication by iw in 
Fonrier space (with i being the imaginary nnit). The signals of the magnetic probes were 
recorded nsing the shotfile system, exactly the same way as dnring a plasma discharge. 
The measnrement of the inpnt and ontpnt signals was started at random phases dnring 
the sweep, bnt the 10 second recording length ensnred that the whole freqnency range is 
covered. The raw signal of one of the cnrrent clamp measnrements is shown in fignre 6(a). 
Figure 6(b) shows the spectrogram for a part of the signal. In this case the sweep started 
at about f = 1.6 sec and the 4 kHz frequency step is clearly visible on the spectrogram. 
Similarly, the raw signal of the magnetic probe and the corresponding spectrogram is 
shown in figures 6(c-d). 
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Figure 6. Raw signals and spectrograms of the signals of transfer function 
measurements, (a) The raw signal of the current clamp measurement of the B31- 
14 coil, (b) The spectrogram of a part of the current clamp measurement, (c) The 
raw signal of the B31-14 coil, (d) The spectrogram of a part of the magnetic probe 
signal. 


The length of a constant frequency period is about 36 ms. For the purpose of the 
present evaluation a 20 ms long interval was selected from the middle of the constant 
frequency periods to safely exclude the frequency jumps. Each 20 ms long interval 
was handled as a separate measurement at a given frequency. The exact frequency of 
the driving signal during the constant frequency periods was measured by finding the 
spectral maximum of each 20 ms long signals. 
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If the changes of the magnetic held at the coil would be known, the ideal transfer 
function could be calculated in the following way: 

^ iuCdBicj) ’ 

where B{u) denotes the Fourier transform of the magnetic held at the coil, 0{u) denotes 
the Fourier transform of the output signal registered by the ADC of the magnetic coil. 
However, B{u) is not know, only the signal of the current clamp is registered which is 
denoted by I{u) since this signal is considered as the input of the measurement. The 
relation between the magnetic held at the coil {B{u)) and the signal of the current clamp 
(/(a;)) is given by the transfer function of the current clamp measurement system: 


= 


B{u) 

Since is not known. 




0{uj) 


( 12 ) 


(13) 


function was calculated separately for the constant Uj frequency periods which, when 
assembled, gave a hrst approximation for transfer function of a given magnetic coil at 
the discrete ojj frequencies. 

The absolute value and phase of (dehned in (13)) of coil B31-14 is shown 

in hgure 7. The absolute value changes more than four orders of magnitudes over the 
0 — 1 MHz frequency range. It has a peak between 20 kHz and 40 kHz which arises from 
the fact that despite of the high amplitude on the magnetic coil, the signal of the current 
clamp strongly decreased in this frequency range. This phenomenon is visible both on 
the raw signals and on the spectrograms in hgure 6 where the relevant time interval is 
indicated with dashed blue lines. Such behaviour of the magnetic coils has not been 
recognized in plasma discharges which suggests that the driving magnetic held strength 
was not measured accurately in this frequency range. This low frequency resonance is 
not the feature of the magnetic coil, but it is most probably due to a resonance in the 
current clamp measurement set-up. Since (dehned in (12)) is not known, but it 

has a signihcant impact on the results, the absolute calibration of the coils is not possible 
from the recorded data. However, the response of the diherent coils was measured by 
using the same setting, thus the measurement is excellent for relative calibration which 
is shown in the next part of this section. 

In order to prove the pertinence of the relative calibration, the relation between the 
ideal transfer function of the system (dehned in (11)) and the transfer function dehned 
in (13) has to be examined. Using (11) and (12) we dehne the relative transfer function 
between two diherent magnetic coils in the following way: 

{u) Ofc(a;)/ia;cdHfc(a;) 






Hf{uj) Oi{u)/iujCdBi{u) 


( 14 ) 
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Figure 7. The absolute value and phase of the calculated transfer function of the 
B31-14 coil. 


Since the same set-up was used for each measurement, the transfer function of the 
current clamp is the same everywhere: H’jfiu) = Hfiu) = Thus the relative 

transfer function can be determined without the knowledge of 




Ok{u:)/Ik{u:) 


(15) 


Oi{uj)/Ii{uj) 

The measurement has been performed on each ballooning coil and their relative 
transfer functions defined in (15) between all ballooning coil pairs were calculated. The 
phases of the relative transfer functions of independent pairs relative to B31-14 are 
shown in hgure 8. The effect caused by the differences in the individual pick-up coil 
transfer functions is most signihcant in the 50-250 kHz frequency interval, which is the 
range of ICRH ion-driven TAE modes and other fast particle driven Alfvenic instabilities 
on AUG [22]. In this interval, the relative phase delay between the different probes can 
reach 1 radian (~ 60°). The signihcant variation in the transfer functions between the 
various probes is most probably caused by effects of the induced mirror currents within 
the in-vessel structures, which is not easily accessible theoretically and by the different 
resistance of the probes (from 57.9 G to 191.4 G). 


5. Correction for the coil response in mode number analysis 

The relative phase differences presented in section 3 can be corrected in the mode 
number analysis by using the relative transfer function dehned in (15). This correction 
signihcantly reduces the systematic errors recognized in the phase of the magnetic probe 
signals. 

For a pure sinusoidal mode structure, the relative phase of the mode described by 
(7) is equivalent to the short time Fourier cross-transform (dehned in (3)) calculated 
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Frequency [kHz] 

Figure 8. The relative phase difference between the transfers functions of ballooning 
coils relative to B31-14. 


from the magnetic field fluctuations: 

ipki = SiTg{Sbk{t,u)Sb*{t,u)} , (16) 

where Sb{t, u) is the short time Fourier transform of the magnetic field in the centre of 
one coil. If Sbi{t,uj) is non-zero - which is a necessary condition - then: 

_ J sb,{t,u) \ 

^[Sbi{t,u)j ■ ^ ’ 

Recalling (11) and considering (14) and (15) ipki can be expressed with the So(t,oj) 
short time Fourier transforms of the measured signals and the Hki{u}) relative transfer 
function: 

J Sok{t,uj)/iuHf{uj) 

( Soi[t,uj)/iu:Hf[ijj) 

which means that applying the correction of the relative transfer functions of the 
magnetic probes evaluated in section 4, the ipki cross-phase of the magnetic fluctuations 
can be more accurately determined from the measured signals. 

The toroidal mode number of the TAE presented in section 3 is evaluated again 
using (18). The resulting toroidal mode number of n = —4 agrees with the result 
achieved without the correction. However, the correction leads to a signihcantly better 
£t, as illustrated in hgure 9a and quantitatively characterized by the differences between 
the residuals shown in figures 3b and 9b. 

The deviation of the relative phases from the n = —4 for independent channel pairs 
as a function of time is visible in figure 10. Contrary to the results shown in figure 3b 
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Toroidal mode number fit Residuals - 1.32 s (g) 144.2 kHz 



Figure 9. (a) The relative phases between all pairs of signals are plotted as a function 
of relative probe position. The different colours represent the different time points as 
shown on figure 2. The best fitting line with n = —4 slope is plotted with black 
solid lines. Compared to figure 3, here the cross-phases were corrected with the values 
derived from the coil response measurement, (b) The residuals (defined in (8)), which 
quantify the differences between the fit and the measurements, are shown as a function 
of the mode number. Compared to figure 3b the n = —4 mode number fits very much 
better. 


and figure 4, here n = —4 is the best htting value even in the [—100,100] interval. 
Without correction, the residual of the best £t (n = —4) in the interval [—20, 20] is 
only 2-3 times lower than the residual of close integer numbers such as —10, —16 and 
8. Using the correction, the residual of the best fit is lower by more than an order of 
magnitude than the second best one. 

The same comparison with and without cross-phase correction has been performed 
for different AUG discharges. In total 11 cases are presented from the period 2006-2013. 
The results are presented in table 1. In most cases an order of magnitude improvement 
is observed in the residual (dehned in (8)) of the fits. However, the remaining deviation 
still has a systematic component. Considering the 11 cases presented in this paper, 
the remaining systematic error has no clear, direct dependence on the frequency or on 
the toroidal mode number. We have to note that on AUG there is a slight difference 
(~ 0.01 rad) in the poloidal position of the toroidal array of balloning coils which - 
taking the poloidal mode numbers into consideration - contributes to the remaining 
systematic error. Some systematic errors may arise due to the angular orientation of 
the pick-up coils or due to unknown inaccuracy in the coil position. Eddy currents in 
the vessel structures can cause major systematic errors. Their exact behaviour could 
only be taken into consideration by simulations [6]. 

An important corollary of the systematic error elimination is that one can 
investigate mode numbers with fewer detectors. This has been demonstrated by choosing 
two coils with a relative probe position ^/tt = 0.131, so that ^/vr is not a low order 
rational. The residual of the linear £t as a function of mode number is shown in figure 11. 
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Figure 10. The deviation of the relative phases from the n = —4 mode for independent 
channel pairs as a function of time. The cross-phase correction leads to a significantly 
better fit compared to the one of figure 4. 
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Time range 

Frequency 

Mode 

Average Q w/o 

Average Q with 

number 

[s] 

range [kHz] 

number 

correction [rad^] 

correction [rad^] 

#21030 

1.68-1.94 

160-180 

-4 

5.39 

0.36 

#21030 

1.79-1.94 

170-190 

-5 

4.53 

0.33 

#23824 

1.27-1.39 

130-160 

-4 

7.91 

0.29 

#23824 

1.22-1.30 

115-135 

-3 

9.17 

0.97 

#23824 

1.15-1.19 

115-125 

-4 

8.62 

0.23 

#23824 

1.30-1.45 

140-180 

-5 

6.46 

0.39 

#26611 

1.60-1.73 

160-200 

-3 

3.97 

1.44 

#26611 

1.60-1.73 

170-210 

-4 

10.71 

7.83 

#26611 

1.63-1.73 

180-210 

-5 

4.39 

0.65 

#29586 

1.16-1.24 

180-220 

-3 

2.27 

0.51 

#29586 

1.16-1.24 

230-190 

-4 

1.33 

0.83 


Table 1. The average residuals (Q) defined in (8) of mode number fits for the 
investigated shots. In most cases an order of magnitude improvement is achieved by 
the correction of cross-phases. 


It shows that even two well positioned probes can provide accurate mode numbers if 
the magnetic held huctuation measurement is precise. However, since (f/n is close to 
2/15 ~ 0.133 the residual of hts with mode numbers smaller or greater than 4 by the 
multiples of 15 are also small as is visible in hgure 11. 

The use of fewer probes can be benehcial in cases when the number of available 
coils is limited. This could happen, for example, when some of the signals are saturated. 
The spacing of the two coils affects the range of mode numbers which can be observed. 
As expressed by Hole et. ah [23], large coil spacing provides hne resolution, but in order 
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Residuals - 1.36 s @ 150.3 kHz 



Mode number 

Figure 11. The residuals defined in (8) are shown as a function of the mode number 
when only two probes (namely B31-01 and B31-02) were used to characterize the 
mode structure. With accurate measurements, fewer detectors are sufficient to provide 
adequate mode number determination. The periodicity of the residuals is because the 
probe placing has a near fifteen fold symmetry {f/n = 0.131 ~ 2/15). 

to identify high mode numbers small coil spacing is needed. It is important to note 
that a spatial aliasing could arise due to the symmetries in the coil spacing. Similarly 
to the case of digitizing time signals, a kind of analogue hlter is needed to eliminate 
the aliasing effect. A cut-off arising from the coil dimensions can serve as an analogue 
hlter, since it cannot detect modes with smaller wavelength than its size. For a high, 
but not unreasonable toroidal mode number {n 50), the coil dimensions would be 
unacceptable. Thus, in a complex set-up with several probes, at least two coils with 
small spacing are suggested in order to decide whether mode number is real or an artefact 
of spatial aliasing. 

6. Conclusion 

In this paper we describe the effect of the magnetic coil transfer functions on mode 
number analysis. The correct characterization of the coil phase transfer functions 
requires in-situ measurements, since this is the only straightforward way to account 
for all the effects in the measurement set-up used in the experiments. A new evaluation 
of previously recorded in-situ coil response measurement data [14] on ASDEX Upgrade 
has shown that the transfer functions of the magnetic pick-up coils differ signihcantly in 
the Alfvenic frequency range. Determining the absolute phase shift between the driving 
perturbation and the pick-up coils is difficult, but for mode number analysis only the 
relative phase differences are important. The phases of the relative transfer functions 
were calculated, and we found that the difference in the relative phase delay can reach up 
to 1 radian (60°) in the 50-250 kHz frequency range. It was demonstrated that applying 
the correction of the transfer function in the mode number analysis leads to considerably 
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better mode number fitting. With this reduction of systematic errors, mode numbers 
can be accurately determined using fewer detectors. The correction method was tested 
on the toroidal mode number of strong, globally coherent Toroidicity-induced Alfven 
Eigenmodes. Eleven example cases are presented from the period 2006-2013. In most 
cases an order of magnitude improvement is observed in the residual of the mode number 
hts. 
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